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The error probability analysis of a noncoherent Differential Phase Shift Keyed 
(DPSK) receiver employing diversity combining techniques is performed. It is assumed 
that the system operates over a frequency-non-selective, slowly-fading Rician channel. 

This thesis analyzes Equal Gain Combining (EGC), Selection Combining (SC) 
and Post Detection Selection Combining (PDSC). The first two diversity combining 
techniques are widely used in communication systems, while Post Detection Selection 
Combining is a new technique. Previous analysis of the EGC and the SC techniques 
shows that the EGC technique has a better performance than the SC technique in a 
Rayleigh fading channel. In this thesis, the effect on the performance of a noncoherent 
DPSK receiver using the diversity combining techniques for Rician fading is examined. It 
is shown that the PDSC technique provides a performance that is better than the SC but 
worse than the EGC technique. The PDSC technique allows a relatively simple receiver 
structure independent of the number of diversity branches. 
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I. INTRODUCTION 



To improve the performance in communication systems, especially in a fading or 
jamming environment, diversity has long been used as an effective technique by reducing 
the effects of channel fading. Fading is a destructive addition of the signal components, 
arriving at the receiver by different propagation paths. When the communications signal 
is subject to fading, the error probability increases due to the amplitude variations in the 
received signal components which are the result of phase variations. The phase variation 
is caused by the time- varying multipath characteristics of the channel. The phase 
differences in the received signal components cause a noncoherent combining loss at the 
receiver. To minimize performance losses one could require to increase the transmitter 
power but in most cases this is undesirable. To reduce the degradation effects of channel 
fading, one of several well-known diversity techniques can be applied either at the 
transmitter or at the receiver. 

Channel diversity is a powerful technique that provides relatively low cost link 
improvement. The basic concept is to provide L independent replicas of the transmitted 
signal over independent possibly fading channels to the receiver. The probability that all 
signals will fade simultaneously is small. One of the transmission paths may be subject to 
a deep fade while some other paths may have strong signal. In such an environment, by 
using more than one signal, both the average and the instantaneous signal-to-noise (SNR) 
at the receiver can be improved. Time diversity of order L is the redundancy where each 
information bit is transmitted L times or once over L channels. Each one of the L 
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channels is a so called a diversity branch which in our case is subject to independent 
Rician fading. 

In general, the diversity methods can be classified as frequency, time, space 
(antenna), angle or polarization diversity. When using frequency diversity, L different 
frequencies carry identical signals. The separation between the carriers must be 
sufficiently large to allow a separation that equals or exceeds the coherence bandwidth 
(A/) c of the channel. Here coherence bandwidth is defined as the frequency range over 

which the signal can pass without distortion. A frequency non-selective channel provides 
a coherence bandwidth that is larger than the signal bandwidth. Another diversity method 
is the time diversity in which L information-bearing signals are transmitted in L different 
time slots. The separation between successive time slots must equal or exceed the 
coherence time (A/) c of the channel. Here coherence time is defined as 

a) 

B d 

where B d is the Doppler spread of the channel. A slowly fading channel provides a larger 
coherence time. 

Space diversity (antenna diversity), which is widely used because of its simple and 
economic implementation, uses a single transmitting antenna and several receiving 
antennas. The receiving antennas must be placed sufficiently apart to have significantly 
different propagation delays to each antenna. Generally a 10-wavelength separation 
between any two antennas is considered to be sufficient. Angle diversity (direction 
diversity) requires a set of directive antennas where each one responds independently to a 
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wave to produce an uncorrelated faded signal. In Polarization diversity, only two diversity 
branches are available. This is achieved when signals that are transmitted use two 
orthogonal polarizations. 

The advantage of frequency and time diversity compared with space, angle and 
polarization diversity techniques is that the number of transmitting and receiving antennas 
does not need to be increased. On the other hand, they require a wider bandwidth. 

To obtain a performance improvement, the information from the diversity 
branches must be combined efficiently. For this purpose a number of powerful and robust 
diversity combining methods is examined to combine uncorrelated fading signals 
obtained from the diversity branches. 

Diversity combining can be accomplished before or after the signal detection. If 
the diversity combining is applied before the signal detection, it is called pre-detection 
combining. Conversely, if it is done after the signal detection, it is called post-detection 
combining. 

In reference 1, equal gain combining and different orders of selection combining 
techniques (SC, SC-2 and SC-3) are investigated for noncoherent Differential Phase Shift 
Keying (DPSK) in a Rayleigh fading channel. The bit-error rate performance of the 
diversity combining methods shows that combining the two or three largest signals (SC-2 
or SC-3) offers significant performance improvement over just selecting the largest signal 
component and the performance improves as L increases. 
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In this thesis, Equal Gain Combining (EGC), Selection Combining (SC) and Post 
Detection Selection Combining (PDSC) techniques are analyzed for the noncoherent 



DPSK receiver [2] in a Rician fading channel. 

The communication channel is modeled as a frequency-non-selective, slowly- 
fading Rician channel and is assumed to be perturbed by additive white Gaussian noise. 

A. RICIAN FADING CHANNEL 

For the receiver investigated in this thesis, the received signal amplitude is 
assumed to be a Rician random variable. As the signal travels to the receiver, it is 
reflected by the surrounding objects which causes indirect (diffuse) signal components. 
The receiver may pick up the direct signal component from the line-of-sight path and an 
indirect (diffuse) signal component from a reflected path. The indirect signal component 
is also referred to as multipath. The multipath components destructively interact and 
reduce the performance of the receiver. This phenomenon is called multipath fading. 

If the line-of-sight component is not zero, then the received signal amplitude is 
modeled as a Rician random variable and the channel is referred to as a Rician fading 
channel. 

For a Rician faded signal, the probability density function of the received signal 
amplitude is [3] 



r 





> * 



( 2 ) 



0 



, otherwise 
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where a 2 is the average power of the unfaded (direct) signal component, 2 cr is the 
average power of the faded (diffuse) signal component and I 0 (») is the zeroth-order 

modified Bessel function of the first kind. The total power of the received signal is the 

sum of the power of the direct and diffuse signal components, a 2 + 2<r 2 . The ratio of the 

direct signal component power to the diffuse signal component power is 

_ CC 2 _ direct component power 
2<j 2 diffuse component power 

Here y characterizes the strength of the fading channel. We assume that each diversity 

reception fades independently and the channel can be modelled as a slowly fading, 
frequency non-selective Rician fading channel. 

When the direct signal component is zero, that is a = 0 and 7 = 0, the received 
signal amplitude is modeled as a Rayleigh random variable and the channel is called a 
Rayleigh fading channel. The probability density function of the received signal 
amplitude in a Rayleigh fading channel is [3] 



M r ) = 



“T e xp(--— 7 ) , r >0 



O' 

0 



20 i 



, otherwise 



(4) 



B. NONCOHERENT DPSK RECEIVER 

The block diagram of the noncoherent DPSK receiver, as used in this thesis, is 
shown in Fig. 1 [2]. 
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Figure 1. Noncoherent DPSK receiver 




DPSK is an efficient modulation method in fading channels when it is difficult 
tosynchronize the phase. It requires phase stability over only two consecutive signaling 



intervals. 

In DPSK, when the data bit £> =0 is transmitted, the waveform v^(f) for two 
consecutive differentially encoded bits c i _ l c i is [2] 

v (1) (r) = ±A p T (t-(i-l)T) cos (2^ f c t+6 0 )± A p T ( t-iT ) cos(2^jf c ?+0 o ) 



= [± A p T (t—(i—l)T)±A p T (t-iT)] cos (27if c t + 6 0 ), 
where p T (*)is a signal pulse of duration of one bit interval given by 

fl , iT<t<(i + l)T)} 



(5) 



p T (t-iT)=- 



0 , otherwise 



( 6 ) 



0 0 is the signal phase, f c is the carrier frequency, A is the amplitude and ± denotes the 
polarity of the two consecutive bits. 

When the data bit b t = 1 is transmitted, there is a phase change of K radians 
between two consecutive bits. Hence, the waveform v^(t) is 

v (2) (f) = ± A p T (t-(i-i)T) cos (2 K f c t+6 0 )+ Ap T (t-iT) cos (2 7tf c t+6 0 ) 



= [± Ap T (t-(i-\)T)+A p T (t-iT)]cos(27tf c t+e 0 ). 



(7) 



If we assume b x = 0 is transmitted, the received signal with white Gaussian noise n(t) is 



of the form 



r(t) = [±A p T (t) ±Ap T ( t-T )] cos (27tf c t +d 0 )+n(t). 



( 8 ) 



The in phase outputs of the receiver are given by 
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(9) 



Y u = 



r at . . , o 

±—cosd 0 +N lc 


+ 


(at 

±cos9+N 


2 ) 




2 ) 



=±ATcos6 0 +N lc 



and 



f AT 

r*= ±— cose 0 +w lt 






( AT 



cos 6 0 +N lc 



= N 



2 c 5 



( 10 ) 



where 



N u =N ic '+Nj and N 2c = N lc ' - N lc " . 



(ID 



/ // 

N lc and N lc are described as 



N 



lc 



T 

J n(t)cos(27rf c t)dt 



o 



and 



N lc = j n(t)cos(2jtf c t)dt . 

T 



( 12 ) 



(13) 



The quadrature phase outputs are given by 



Yu = 



AT ' \ 

±— -(-sin d 0 ) + N u + 



AT 

± (-sin0 o ) + W, 



=+AT sin 6 0 + N u (14) 



and 



Y 2s = 



AT ' ^ 

±— (-sin e 0 ) + N ls 



( AT 

±— (-sin 0 O )+W„ 



= N 



2s ’ 



where 



N ls =N u +N ls and N 2s =N u -N u . 



(15) 



(16) 



/ // 

N [s and N u are described as 
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(17) 



/ 



N 



1j 



T 

= Jn(? ) sin (2^C c r )dt 

o 



and 

2 T 

N u = jn(t)sin{27tf c t)dt. (18) 

T 

The noise n(t) is modeled as additive white Gaussian noise with a two sided power 

N 

spectral density . The noise is assumed to be unaffected by the fading channel. 
Therefore N lc , N 2c , N ls and N 2s are independent, identically distributed (i.i.d.), zero 



mean Gaussian random variables with variance 



NJ_ 

2 



If we assume that b , = 0 is transmitted, then for the lower branch (Fig.l.) which is the 
signal branch, the outputs are given by 



Y u =±ATcosd + N lc 
Y u = +AT sin 6-+N ls . 



(19) 



The upper branch, which is the non-signal branch, is described by 

Y =N 

l 2c ■' v 2c 

Y =N 

1 2s •' V 2i ' 

The decision variables are defined as 

fl, 2 =Y u 2 +Y u 2 ={±ATcosd + N ic f +{+ATsmd + N u ) 2 
for the signal branch and 

R 2 2 = y 2c 2 + y 2s 2 = N2 C 2 +n 2s 2 

for the non-signal branch. 



( 20 ) 



( 21 ) 



( 22 ) 
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In chapter H, HI and IV the bit error probability expressions of the EGC, the SC 
and the PDSC techniques are analytically obtained and the numerical results are 
presented in chapter V. 
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II. EQUAL GAIN COMBINING 



Equal Gain Combining (EGC) is a technique in which the information from all 

available diversity branches (paths) is added with equal weights. In other words, the L- 

branch equal gain combiner outputs V, which is the decision variable for the signal 

branch, and V 2 which is the decision variable for the non-signal branch, are the sum of L 

independent random variables, 

v t =R n 2 + R l2 2 +R l3 2 +... + R lL 2 
v 2 =R 2 \ 2 +R -22 +^ 23 2 + -~ + R 2 L » 

where R u 2 and R 2i 2 , i = (1,2,.. .L) , are the amplitudes of the signal received by the i' h 

channel. The receiver for EGC is shown in Fig. 2. 

This method has been commonly used by receivers that use noncoherent 
demodulation and is considered to be equivalent to Maximum Ratio Combining (MRC) 
for coherent demodulation [3]. 

Due to the “noncoherent combining loss” effect, the EGC method is not optimal. 
In particular, increasing the number of diversity channels (L) does not necessarily 
improve the performance. 

A. PROBABILITY DENSITY FUNCTION OF THE DECISION VARIABLES 

In the noncoherent DPSK receiver as shown in Fig. 2, the decision variable V, is 
the sum of the squares of 2L Gaussian random variables. If we assume that the data bit 
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Figure 2. Noncoherent DPSK receiver for Equal Gain Combining (EGC) with L‘ h order 
diversity 
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b i = 0 is transmitted, then the output for the signal present branch is 




( 24 ) 



For convenience, we normalize (21) as follows; 









-|2 



VF 



Jj 






sinS + ^ 



->2 



4t 



L 

2 



± -v/2 T A cosd 



VF 



VF 



J 



(25) 



/ 



2 [2” [2" ^2r 

By setting i? u 2 =-i?,\ n lc< = ^|— W le , n 1% = ^-(V u , A, = ± — A and 0, =0, 



we 



get 

-^i* 2 = (2A t cos 9 k +n XCt ) 2 +(— 2A t sin 0*+« 1Jt ) 2 . (26) 

For a Rician fading channel A k is a Rician random variable with parameters a and o. 
A k cos 0 k is a Gaussian random variable with mean a c -a cos <p where <p is an arbitrary 
phase and variance a 2 , A k sin Q k is a Gaussian random variable with mean 
a s — —a sin (j) and variance <7 2 , 2 A k cos 6 k is a Gaussian random variable with mean 
cc c '= 2 a cos 4> and variance a' 2 = 4a 2 and 2 A k sin 6 k is a Gaussian random variable with 
mean a s ' = -2asin <p and variance o' 2 = 4cr 2 [4]. All of these L independent random 
variables 2 A k cosd k for k=l, 2,..., L have a mean of a c ' and all of the L independent 
random variables 2A k smd k for k=l, 2,..., L have a mean of a/. The zero mean 
Gaussian random variables n lc t and n Uk for k=l, 2,..., L are independent with variance 



13 



o n 2 = N 0 . Therefore the probability density function of V, is non-central Chi-square 
distribution with 2L degrees of freedom [3], 



1 



/V ' (V,) = 2 
2<Jj 



L-\ 

v i Y*"" 

— exp 


' (v,+s 2 )^ 


1 L-\ 


2 


[ s ) 


2a, ) 




J 



, v, >0 , (27) 



where 



a 2 =o' 2 +a n 2 — 4(7 2 +<7 2 



(28) 



and 



£. z. 

5 2 =X« C ’ 2 +X a / 2 = L(2a cos (f>) 2 + L(- 2a sin <f>) 2 =4La 2 . (29) 

i=i i=i 



We denote the average energy per diversity channel as 

E = a 2 +2<j 2 . 

The bit energy is related to the energy per diversity channel by 

X = LE => £ = — . 



(30) 



(31) 



For the non-signal branch, the decision variable V 2 is sum of the squares of 2L Gaussian 
random variables with zero mean and variance o 2 . Hence, the non-signal branch output 
is given by 



V 2=^Z R 2 * 2 = X i n 2c t 2 + n 2s t 2 )> 



(32) 



where 



k= 1 k = I 



/? 2* 2 =«2c t 2+ «2, t 2 



(33) 



and 
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CJ 2 2 = var(/i 2ct 2 ) = var(« 2jt 2 ) = a n 2 = N 0 . 



(34) 



Hence the probability density function of V 2 is Chi-Square with 2L degrees of freedom 
which is given by [3] 



/v,( V 2 ) = 



(2<t ! 2 )‘(L-1)! 



L - 1 



exp 



2(7, 



v 2 >0 



(35) 



B. BIT ERROR PROBABILITY 



We now have the decision variables for both the signal and the non-signal 
branches under the assumption that b t = 0 is transmitted. The conditional probability of 
correct detection is defined by [2] 

Pr{ Correct Detection I V, = v, }= Pr (V 2 < V l I V, = v, ) . (36) 

Conversely, when b- t = 1 is transmitted, due to the symmetry of the receiver branches, the 
conditional probability of correct detection is described as 

Pr{ Correct Detection I V 2 = v 2 }= Pr (V, <V 2 \V 2 = v 2 ) . (37) 

Writing out (36) we obtain 



Pr{V 2 <y, I V,=v,}=j/ Vi (»,)*>, 

0 




1 

(2ct 2 2 ) l (L-1)! 




exp 

V 



2(7, 



dv. 



(38) 



By using the identity [4] 
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, u,n> 0 (39) 



— n ! u k 



J*" exp (~fJx)dx = —^ r - exp (-/*«)£ — — 

J 0 jU ^ k\ ji n 



and setting u = v, , n — L-\ and /J. - — , we get 

2a, 



Pr{ V 2 < V. I V. = v, }= 

(2a 2 2 ) L (L-l)! 



X 



q-D! 

1 * 



-exp - 



2a 2 , 
v 2 / 



Z .-1 



2a, 



(L-l)! 






Jk! ( , V"* 



v 2 ^ 2 ; 



= l-exp 



^ L-1 J / .. V 



2c T. 



J 



y- 



2CTn 



(40) 



This result is inserted into the probability of correct detection expression given by 

P c = Pr{ Correct Detection }= |Pr{ V 2 <V, I V, =v, }/ Vi (v,) dv x (41) 



=1 



1-exp 



f ^ 

v, 



2c\ 



L - 1 1 f 

Y- 
6 *! 



\* 



2 ( 7 , 



2 ( 7 , 



(v,+*’)) 


J 


rv^i 


V J 


1 L-l 


2 

<J X 

\ 1 / 



xexp 



Expanding this expression leads to 



dv, . 



? 1 


L-l 

(v. XT 


(v,^ 2 )l 






J 0 2a, 2 


2 eX P 

l 5 J 


2a, 2 J 


l -, 2 J 



dv, 



oo 

-Jexp 



2 ( 7 , 



''L-l 1 f 

Z- 



V 



i-l 



2a, , 
2 / 



2 ( 7 / 



fv, Yi 



v* 2 , 



exp 



L-l 

(„ XT 



\° / 



(42) 






(43) 



' ^ n *, 2 



dv,. 
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In order to simplify, we use the identities [5] 

i) I n (z)={-jT J n Uz\ 

ii) ]x m+2 cxp{-ax)J n {2p^fx)dx = -^^ r exp 



( 



i«) K 



P 



a ; 



m 

-X 



to: 



fm + n ^ 



y m- P) 



/? 2 Y 





T n 




< a j 


L m 


< 01 J 






(44) 

(45) 

(46) 



^5 p: 

The first integral expression in (43) is the integral of a probability density function given 
by 



u 



L-\ , 

'v, ^ 



\ s ) 



0 2(7, 

The second integral becomes 



exp 



2(7, 



h+s 2 )) 

~ 2 1 1 - 1 2 



dv x = 1 . 



(47) 



\ 1 



00 ( 

Jexp - 



V, 



\ 



2(7 o 



L-l 1 

E 1 

0 



.. Y 



2(7, 



2(7/ 



^ rr ( (v l+ s>) 



exp 






^ / — > \ 
V v i 5 



1 



^ 2 v- 

,+ ^ 

G 2 



vC-O" 



exp 



2(7/ 



1- 



2 . 2 

(7, + (7 2 



2(7, 

2 v\ w r 

X 

;J k ~° \ 



dv x = 



i + ^ 



2 \ 



-k 



(48) 



m=0 



ml 



pfc + L-1^ 


-s 2 (7 2 


k-m J 


2g 2 (7, 2 +(7 2 2 



Substituting (47) and (48) into (43), the probability of correct detection is obtained as 
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p= 1 - 
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^ 2 
1 + — 



r ex p 



2C7, 



2 

1- a2 



^ L-! f 



2 2 

CT, +(7 2 



> ) 
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1 + — 



2 "\ 



-k 
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( 49 ) 



The probability of error expression is obtained from P e = 1 — P c . So 



P =1-P = — 

e x A c , 



1 + — 
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(50) 



From (28), (29) and (34), we can express <7, 2 , s 2 and a 2 as follows; a / =4<t 2 + /v r 0 , 
,? 2 =4La 2 and a 2 2 =N 0 , respectively. Then using the average energy per diversity 
channel E in (30), we obtain 



E = a z +2a* => 2cr 2 = 



2 E 



7+1 



(51) 



where y-a 2 12a 1 is the direct-to-diffuse ratio given in (3). The probability of error in 
terms of the average energy per diversity channel is given by 
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P e = exp 
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By using (3 1 ), the probability of bit error for a noncoherent DPSK receiver with equal 
gain combining in a Rician fading channel in terms of the bit energy-to-noise density ratio 



E b / N 0 is given by 



P b =exp 
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(53) 



To verify this result, we can set L=1 and s 2 = 4La 2 = 0 (i.e.; there is no direct signal 
component) and we obtain the well-known expression of the probability of bit error for 
DPSK in Rayleigh fading channel which is given by [3] 



Pb 2 E b 
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( 54 ) 
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Similarly, if we set L=l, we obtain the well known expression for probability of bit error 
of DPSK in Rician fading channel given by [6] 
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numerical results of this analysis are presented in 
SC technique is analyzed and the bit error probability 



chapter V. 
expression 



( 55 ) 



In the next 
is obtained. 
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III. SELECTION COMBINING 



Selection Combining (SC) is a predetection combining technique that takes the 
signal that has the largest amplitude (hence largest signal-to-noise ratio) in the L branches 
at the output of the combiner preceeding the DPSK demodulator. We define 

V = max(Y 1 ,Y 2 ,Y i ,...,Y L ), (56) 

where Y k , (k = l,2,...,L) , is the squared amplitude of the signal received by the k ,h 

channel. The receiver for SC is shown in Fig. 3. The derivation of the probability density 
function for the largest random variable V is given in the Appendix. 



A. PROBABILITY DENSITY FUNCTION OF THE DECISION VARIABLES 

The pdf of Y k is non-central Chi-square and is given by [3] 



2cT 



G , * + * 2 )l 


J 
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l o 


G 2 

\ / 



; 



The cumulative distribution function (cdf) of this pdf is [3] 



F r .(y,)=l-Q 



_2 ’ 2 

G G 



V 

Q(a,b) in (58) is known as Marcum’s Q function and defined as [3] 

f (x 2 + a 2 )^ 



Q{a,b)= Jxexp 



I 0 (ax)dx 



J 



(57) 



(58) 



(59) 
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Figure 3. Selection Combining with L' h order diversity 



and 



Q{a,0)= 1 



Q{0,b)=e\p 



b 2 ^ 



(60) 



Now we can define the probability density function for the largest diversity branch by 
using (A. 14) as derived in the Appendix [7] 

f v {v)=L f Yk (v)F yt L ~ x ( v ) . (61) 



Substituting (57) and (58) into (61), we get 



2cT 



r (v + S 2 )^ 




f r \ 

Vv s 


1 — (1 


s £) 


2<r 2 J 


h 


a 2 

\ / 


1 


,°’ a 1 



L - 1 



(62) 



B. BIT ERROR PROBABILITY 

In deriving the bit-error performance, the bit error rate expression of DPSK for a 
non-fading channel is used which is given by [3] 
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_v N 
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(63) 



where v/N 0 is the signal-to-noise ratio per channel at the input of the DPSK 

demodulator. To obtain the error probability expression, we evaluate an integral of the 
form 



p>=]p t w/vM to- 
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Substituting (62) and (63) into (64), the probability of bit error expression for SC leads to 
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Changing the variable, u = — in (65), we get 
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The direct-to-diffuse signal power ratio of the signal is defined as 



( 66 ) 
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(67) 



7 = I+' 



Hence the energy per diversity channel is described by (eq. 30) 
_ . Y c 2 ^ 



E = s + 2a = 2<J‘ 



2a 2 



-+1 



= 2a 2 (y+l) => a 2 = 



2 ( 7 + 1 ) 



( 68 ) 



Rewriting the parameters a 2 and s 2 in (66) by using (67) and (68), the probability of the 
bit error expression in terms of the energy per diversity channel is obtained as 
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(69) 



As a result, the probability of bit error for noncoherent DPSK receiver with Selection 
Combining in a Rician fading channel is defined as 
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To evaluate this expression numerical analysis is required. Some approximations for 
computing Marcum’s Q function can be found in the literature. In this thesis, a recursive 
method [8], is used because of its more stable and accurate results. For the numerical 
evaluation of I 0 (x), the exact definition 



\ 2 ' 
x ' 



/oM=I 



s G'O 2 



(71) 



can be used [9]. 

The numerical results of this analysis are presented in chapter V. In the next 
chapter, the PDSC technique is analyzed and the bit error probability expression is 
obtained. 
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IV. POST DETECTION SELECTION COMBINING 



Post Detection Selection Combining (PDSC) invented by Professors Tri Ha and 
Ralph Hippenstiel [10] is a technique which chooses the largest values from the signal 
branch and the non-signal branch as the decision variables in a DPSK noncoherent 
demodulator as shown in Fig. 4. Between the L diversity branches, the maximum 
amplitudes are chosen independently as the decision variables in both branches. For the 
signal present branch, the decision variable Vj is given by 



For the non-signal branch, which corresponds to the noise only condition, the decision 
variable V 2 is given by 



A. PROBABILITY DENSITY FUNCTION OF THE DECISION VARIABLES 

For the signal branch, the probability density function of the i‘ h branch output F,, 

before combining, is modeled as in (27) with a non-central Chi-square distribution with L 
set to 1 [3]; 





branches. 




( 74 ) 
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Figure 4. Noncoherent DPSK receiver for Post Detection Selection Combining (PDSC) 
with L‘ h order diversity 
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Each of the L diversity branches has the probability density function given by (74) and a 
corresponding cumulative distribution function, cdf, given by 

' VST 






<Ti <7, 



(75) 



Since all L branch variables are statistically independent, using (61), the pdf for the signal 
branch decision variable becomes 



A(v,)=i4(v,)F,/' , (v,) 
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Similarly, in the non-signal branch, under the noise only condition, the maximum 
amplitude is chosen as the decision variable. Before diversity combining, the i ,h branch 
output Y 2i is modeled as in (35) with a Chi-square distribution with L set to 1 : 
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(77) 



with a corresponding cdf given by 

F r„.(>’2,) = l-exp 



2(7 2 2 , 



(78) 



The probability density function for the largest amplitude in the non-signal branch 
becomes 



fv 2 ( V 2 ) = L fy 2i ( V 2 ) F Y 2 , L ~' ( v 2 ) 
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B. BIT ERROR PROBABILITY 

We now have the decision variables for both branches. Using the conditional bit 
error probability expression we get 
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Applying the binomial theorem, 
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to the exponential term in (81), the exponential term can be simplified as 
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and (81) becomes 
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( 84 ) 



L-l f T-\\ 

■L I 



k=0 



(- 0 * t— r ex P 

\ + k 



f 


' 1 +ik 


Y| 


~ V 1 


2 a 2 




\ 


2 


Jj 



Substituting (76) and (84) into (80), we get 
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Rewriting the parameters a 2 , a 2 2 and s 2 as 
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(87) 



and 



s 2 =4a 2 . (88) 

Then the probability of bit error expression in terms of the energy per diversity channel is 
described by 
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By a change of variable, u = — , we get 
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Substituting E by E b / L in (90), we get 
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As in selection combining, this expression is evaluated by numerical analysis. The results 
of this analysis are presented in the next chapter. 
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V. NUMERICAL RESULTS 



In chapter H, HI and IV the bit error rate (BER) expressions for noncoherent 
DPSK signals operating in a frequency non-selective, slowly fading Rician channel for 
EGC, SC and PDSC techniques are analytically obtained. The main objective of this 
thesis is to illustrate the performance by allowing the comparison of the EGC, SC and 
PDSC techniques. The numerical evaluation of the bit error rate expressions is performed 
by using MATLAB 5.1 [11] and is presented in Figures 5-26. The bit energy-to-noise 
ratio is chosen in the range of 6-20 dB. A direct-to-diffuse signal ratio of y = 0 denotes a 
Rayleigh fading channel. A y of 1 is effectively a Rayleigh fading (weak Rician) channel 
while y of 10 is effectively a Rician (i.e.; there is a strong direct signal component) 
channel. Hence, in numerical evaluations, values of y = 0,1,5,10 provide sufficient detail 
to illustrate the performance differences. 

In Figs. 5-10, the EGC, the SC and the PDSC performances of the receiver are 
illustrated separately for diversity orders of L=2, 3, 4, 5. As expected, the L=1 case 
provided the same performance for all three techniques and is not plotted. 

In Figs. 1 1-26, the performance comparison of the EGC, the SC, and the PDSC is 
illustrated for diversity order of L=2, 3, 4, 5 and y = 0, 1, 5, 10 . 

The results in this thesis are based on the assumption that the direct-to-diffuse 
signal ratio y is the same for each diversity branch (i.e.; the average power of the 
diversity branches is assumed to be equal). 
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In Figs. 5-6, the EGC performance of the receiver is shown for diversity orders of 
L=2,3,4,5 and the direct-to-diffuse signal ratios of y = 1,10 . For y = 1 , which means that 
the channel fading is effectively Rayleigh (weak Rician), as L increases the receiver 
performance improves. For y = 10, when the channel fading is effectively Rician (strong 
direct signal component), all BER curves illustrate better performance than y = l case. 
For low values of E b / N 0 , the system performance with smaller diversity order L is 
observed to be superior to those with larger diversity order L. This trend reverses as 
E b / N 0 increases. This effect is due to the noncoherent combining loss. As mentioned 

before, noncoherent systems are not optimal because of this phenomenon. This effect is 
more evident with larger L values as seen in Fig. 6. We observe the performance 
improvement effect of increasing L after a certain bit energy-to-noise density ratio, 
primarily since the noise in each diversity channel is significant at low E b / N 0 values. As 

E b / N 0 increases, the noise in each diversity channel becomes less significant. In the 

EGC technique, as the bit energy-to noise density ratio increases (i.e.; as noise becomes 
less and less significant), each diversity branch contributes positively to the overall SNR 
and therefore a system with a larger diversity order L will perform better than a smaller L. 

Figs. 7 and 9 demonstrate the receiver performance of the SC and the PDSC 
techniques for y = 1 , respectively. As L increases, performance improvement is observed 
due to the weak Rician (i.e.; Rayleigh) fading channel. Since the direct signal component 
is weak, most information is obtained from the diffuse signal components. When a larger 
diversity order L is employed, the probability that the receiver chooses a stronger signal 
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branch becomes larger. In this case, performance improvement as L increases is an 
expected result. For the 7 = 10 case (strong direct signal component), L=2 provides the 
best performance but we observe the slight performance improvement effect by 
increasing L after a certain E b / N 0 value as in the EGC performance for 7 = 10 (Fig. 6 ). 

Again the noncoherent combining loss is observed in both SC and PDSC. Increasing L 
does not provide significant performance improvement since a strong direct signal 
component is already available. 

In Figs. 11-26, the performance comparison of the EGC, the SC and the PDSC 
techniques is illustrated for L=2,3,4,5 and 7 = 0,1,5,10 . As clearly observed in the figures, 
the PDSC technique provides a BER performance better than that of SC and inferior than 
that of EGC. For a particular diversity order L, as 7 increases all three techniques 
demonstrate performance improvement since the receiver performance is greatly 
improved by a strong direct signal component as indicated before. When L is increased, 
performance improvement is observed due to an increasing number of diversity branches 
which permits choosing a higher amplitude diversity branch. In both increasing L and 
increasing 7 cases, the PDSC technique shows a BER performance that is better than the 
SC but worse than the EGC. 
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Bit Error Probability 




Figure 5. The performance of noncoherent DPSK in a Rician fading channel with 
EGC for L=2, 3, 4, 5 and 7 = 1. 
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Figure 6. The performance of noncoherent DPSK in a Rician fading channel with 
EGC for L=2, 3, 4, 5 and 7 = 10 . 
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Figure 7. The performance of noncoherent DPSK in a Rician fading channel with 
SC for L=2, 3, 4, 5 and y = 1 . 
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Figure 8. The performance of noncoherent DPSK in a Rician fading channel with 
SC for L=2, 3, 4, 5 and 7 = 10. 



41 



Bit Error Probability 




Figure 9. The performance of noncoherent DPSK in a Rician fading channel with 
PDSC for L=2, 3, 4, 5 and y = 1 . 
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Figure 10. The performance of noncoherent DPSK in a Rician fading channel with 
PDSC for L=2, 3, 4, 5 and y = 10 . 
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Figure 11. The performance of noncoherent DPSK in a Rician fading channel with 
SC, PDSC and EGC for L=2 and y = 0 . 
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Figure 12. The performance of noncoherent DPSK in a Rician fading channel with 
SC, PDSC and EGC for L=2 and y = 1 . 
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Figure 13. The performance of noncoherent DPSK in a Rician fading channel with 
SC, PDSC and EGC for L=2 and y = 5 . 
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Figure 14. The performance of noncoherent DPSK in a Rician fading channel with 
SC, PDSC and EGC for L=2 and 7 = 10. 
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Figure 15. The performance of noncoherent DPSK in a Rician fading channel with 
SC, PDSC and EGC for L=3 and y = 0 . 
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Figure 16. The performance of noncoherent DPSK in a Rician fading channel with 
SC, PDSC and EGC for L=3 and y = 1 . 
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Figure 17. The performance of noncoherent DPSK in a Rician fading channel with 
SC, PDSC and EGC for L=3 and 7 = 5. 
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Figure 18. The performance of noncoherent DPSK in a Rician fading channel with 
SC, PDSC and EGC for L=3 and 7 = 10. 
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Figure 19. The performance of noncoherent DPSK in a Rician fading channel with 
SC, PDSC and EGC for L=4 and y = 0 . 
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Figure 20. The performance of noncoherent DPSK in a Rician fading channel with 
SC, PDSC and EGC for L=4 and 7 = 1 . 
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Figure 21. The performance of noncoherent DPSK in a Rician fading channel with 
SC, PDSC and EGC for L=4 and y = 5. 
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Figure 22. The performance of noncoherent DPSK in a Rician fading channel with 
SC, PDSC and EGC for L=4 and 7 = 10. 
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Figure 23. The performance of noncoherent DPSK in a Rician fading channel with 
SC, PDSC and EGC for L=5 and y = 0 . 
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Figure 24. The performance of noncoherent DPSK in a Rician fading channel with 
SC, PDSC and EGC for L=5 and y = 1 . 
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Figure 25. The performance of noncoherent DPSK in a Rician fading channel with 
SC, PDSC and EGC for L=5 and y = 5 . 
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Figure 26. The performance of noncoherent DPSK in a Rician fading channel with 
SC, PDSC and EGC for L=5 and 7 = 10. 
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VI. CONCLUSIONS 



In this thesis, the EGC, the SC and the PDSC bit error probability performances of 
a noncoherent DPSK receiver with L' h order diversity have been analyzed in a frequency 
non-selective, slowly fading Rician fading channel. Based on previous analyses of the 
EGC and the SC techniques for a noncoherent DPSK receiver, it is known that the EGC 
technique provides a better performance than the SC method in a Rayleigh fading channel 

[ 13 . 

The EGC technique is widely used in communication systems that use 
noncoherent demodulation, but it has the undesirable features of noncoherent combining 
loss and receiver complexity which depends on the diversity order L. In a Rician fading 
channel and at low values of E b / N 0 , the system performance with small L is superior to 

those with larger L which is due to the noncoherent combining loss. 

The SC technique has been suggested as a simpler technique that can provide 
adequate performance without L dependence. However, it is not an optimal combining 
technique since it does not use all of the available diversity branches simultaneously. 
Similar to EGC, performance improvement is observed when E b / N 0 increases and the 

same effects of noncoherent combining loss are observed. 

Since L may vary with locations as well as time, having receiver complexity 
dependent on L is undesirable. Therefore, it is of interest to find combining techniques 
which are independent of L and which can provide adequate performance. The SC 
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technique is one simple suboptimal design which can achieve acceptable performance 
under typical channel conditions. The EGC and the SC represent the two extremes in 
diversity combining with respect to the number of signals used for demodulation. 

The PDSC technique requires a less complex receiver than the EGC and may be 
implemented regardless of the number of the diversity branches. Furthermore it does not 
require a predetection combiner. As seen from the numerical results, its performance is 
worse than the EGC but it offers a performance improvement over SC. 
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APPENDIX: DERIVATION OF THE PROBABILITY DENSITY 
FUNCTION FOR THE LARGEST RANDOM VARIABLE 



In this appendix, the probability density function f z (z) of the largest random 
variable Z is derived. Lets consider L independent, identically distributed (i.i.d.) random 



variables X V X 2 ,X 3 ,...,X L with the probability density functions 

fx,M = f*i fa ) = /*, fa ) = ••• = fx L ( x l ) = fx M (A. 1 ) 

and we want to select the largest random variable Z 

Z = max{Xj,X 2 ,X 3 ,...X L j. (A.2) 

According to the order statistics, L random samples are arranged in an increasing order as 
X, g)< X „(/;)< X,G)s ... £ X„g)£ ... £ X.fi) . (A.3) 

We next form the L random variables z ; such that 

z, = * rj < z 2 = x f2 < ... < Z k = x rt < ... < z L = x rt . (A.4) 

The k th largest random variable between L random variables is called the k ,h order 
statistics and we determine its density by using 

f k (z)dz-P{z<z k <z + dz} . (A.5) 



The event A = {z< z k < z + occurs if and only if k-1 of the random variables are less 
than z and one is in the interval (z, z+dz). By defining the sets as 

/?!={x<z}, /3 2 ={z< x< z + dz] and /J 3 = {z + dz < *} (A.6) 

the probabilities of events are equal to 
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and 



P(j3 l )=p l =F x {z), 
p (Pi)=P2 =f x {z)dz 



(A.7) 



(A.8) 

p (P3)=P3=l-F x {z + te)- (A.9) 

where F x (x) is the cumulative distribution and f x (x) is the density of x. We can define 
“ p. occurs k i times” as [6] 

— ;Pi*' Pi 1 Pi 1 , +A: 2 + ^3 =L. (A. 10) 

The event A occurs if /?, occurs k-1 times, fi 2 occurs once, and /? 3 occurs L-k times. So 

k x =k- 1, fc 2 =l, k^-L-k (A. 11) 

and by inserting these expressions into (10), we get 

f k (z)dz = P{z<z k <z + dz }= ' ( z )A & ~ ( z + *)] L ’* • 

(A. 12) 

As a result, the order statistics is obtained as 

/»(z) = (*-lKL-*) F “bi 1 - F &T‘ ' f.(z) ■ (A. 13) 

In our case , we want the largest value which means that k=L. So the probability 
density function of the largest random variable between L random variables is defined by 

fz(z)=Lf„(z)F x L -'(z). (A. 14) 
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